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Abstract 

We consider the leptonic asymmetry generation in the vMSM via hadronic decays 
of sterile neutrinos at T <C Tew-* when the masses of two heavier sterile neutrinos are 
between m n and 2 GeV. The choice of upper mass bound is motivated by absence of 
direct experimental searches for singlet fermions with greater mass. We carried out 
computations at zero temperature and ignored the background effects. Combining 
constraints of sufficient value of the leptonic asymmetry for production of dark matter 
particles, condition for sterile neutrino to be out of thermal equilibrium and existing 
experimental data we conclude that it can be satisfied only for mass of heavier sterile 
neutrino in the range 1.4 GeV < M < 2 Gev and only for the case of normal hierarchy 
for active neutrino mass. 

1 Introduction 

The Standard Model (SM) is minimal relativistic field theory, which is able to explain almost 
all particle physics experimental data [1]. However, there are several observable facts, that 
cannot be explained in the SM frame. Firstly, the neutrinos of SM are strictly massless, that 
contradict to the experimental fact of the neutrinos oscillations [2, 3]. The second problem 
is the impossibility to explain the baryon asymmetry of the Universe (BAU) within the SM. 
Finally, the SM does not provide the dark matter (DM) candidate. Also the SM can not 
solve the strong CP problem in particle physics, the primordial perturbations problem and 
the horizon problem in cosmology, etc. 

The solutions of the above mentioned problems of the SM require some new physics 
between the electroweak and the Planck scales. An important challenge for the theoretical 
physics is to see if it is possible to solve them using only the extensions of the SM below the 
electroweak scale [4]. 

The Neutrino Minimal Standard Model (uMSM) is an extension of the SM by three 
massive right-handed neutrinos (sterile neutrinos), which do not take part in the gauge 
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interactions of the SM x . The model was suggested by M. Shaposhnikov and T. Asaka [5, 6]. 
The masses of sterile neutrinos are predicted to be smaller than electroweak scale, and thus 
there is no new energy scale introduced in the theory. The parameters of the vMSM can 
be chosen in order to explain simultaneously the masses of active neutrinos, the nature of 
DM, and BAU. 

The lightest sterile neutrino (the mass is expected to be in the KeV range [1]) can be 
intensively produced in the early Universe and have cosmologically long life-time. So, it 
might be a viable DM candidate. The sufficient amount of this neutrinos can be generated 
through an efficient resonant mechanism proposed by Shi and Fuller [7]. 

In the vMSM the required amount of the leptonic asymmetry (in accordance with Shi 
and Fuller mechanism) can be created due to decays of the two heavier sterile neutrinos. 
This particles are generated at temperature T > T EW and their masses are expected to 
be in range m n < Mj < Tew [8], where m n is the pion mass and Mj is the mass of I- 
sterile neutrino. The leptonic asymmetry at the temperature of the sphaleron freeze-out 
(T ~ Tew) is related to the baryon asymmetry of the Universe. At temperature T < Tew 
the leptonic asymmetry from decays of heavier sterile neutrinos can not convert into the 
baryon asymmetry and is accumulated. As it was shown in [4, 9] the required amount of the 
leptonic asymmetry A = AL/L = {tie — n L )/{riL + ni): 

l(r 3 <A<2/ll (1) 

has to already exist in the Universe at the moment of the beginning of production of the 
DM particles (takes place at the temperature around 0.1 GeV). 

We consider here the leptonic asymmetry generation at T <C Tew, when the masses 
of two heavier sterile neutrinos are between m n and 2 GeV. The motivation is following. 
The mass of heavier sterile neutrino can not be less then m n (the constraint is coming from 
accelerator experiments combined with Big Bang Nucleosynthesis (BBN) bounds [10, 11]) 
and there is no direct experimental searches for singlet fermions with mass more then 2 Gev 
[10]. 

Since the masses of active neutrinos in the vMSM are produced by the "see-saw" mech- 
anism [12] some constraints on the parameters of the vMSM come from active neutrino 
parameters that can be found from the experiments on the neutrino oscillations. Namely, 
these are the mass squared differences of active neutrinos and the mixing angles. Until re- 
cently the mixing angle #13 was supposed to have a close to zero value. But new observations 
indicated its essential difference from zero [13]. 

The aim of this work is to obtain constraints on the parameters of the vMSM from 
the required amount of the leptonic asymmetry and cosmology conditions. Also we want to 
investigate the influence of non-zero mixing angle # 13 on space of the allowed parameters of 
the vMSM. We do it following [14] using a simple model: we ignore the background effects 
and do computations at zero temperature. 

The paper is organized as follows. In Section 2 we present the Lagrangian of the uMSM, 
make its convenient parametrization and present the Yukawa couplings in terms of active 
neutrinos mass matrix parameters. In Section 3 we derive the expression for the leptonic 
asymmetry. The limitations on the vMSM parameters are imposed in Section 4. Section 5 
is devoted to the analysis and conclusions. 

lr This is why these neutrinos are called sterile neutrinos. The left-handed neutrinos of the SM are called 
active neutrinos. 
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2 Basic formalism of the vMSM 



In the vMSM [5, 6] the following terms are added to the Lagrangian of the SM (without 
taking into account the kinetic terms): 

C ad = -F aI L a ®v IR - ^±Lv c m v JR + h.c, (2) 

where index a = e, /i, r corresponds to the active neutrino flavors, indices /, J run from 1 to 
3, L a is for the lepton doublet of the left-handed particles, vm is for the field functions of 
the sterile right-handed neutrinos, the superscript c" means charge conjugation, F a j is for 

the new (neutrino) matrix of the Yukawa constants, Mjj is for the Majorana mass matrix 
of the right-handed neutrinos, $ is for the field of the Higgs doublet, $ = ia 2 ^* . 

After the spontaneous symmetry breaking the field of the Higgs doublet in unitary gauge 

is 



v+h 
V2 

where h is the neutral Higgs field and the parameter v determines minimum of the Higgs 
field potential (v = 247 GeV). In this case Lagrangian (2) acquires the Dirac-Majorana 
neutrino mass terms: 



C DM = —FajUaUm - ^PI r vjr + h.c, (3) 



or in conventional form [15] 

M DM 



C DM = -[ {N L y ^— N L + h.c. ) , (4) 



where 



( M L M J 



M L = 0, M D = F + ^=, M R = M\ (6) 

where M, F are square matrix of the third order with elements F a j and Mjj. 

In zero approximation the vMSM Lagrangian is assumed to be invariant under U(l) e x 
U(l)fj, x U(1) T transformations, that provides preservation of the e, /x, r lepton numbers 
separately. It is also assumed that two heavier sterile neutrinos interact with the active 
neutrinos, but the third (lightest) sterile neutrino does not interact 2 . This assumption can 
be realized by following matrix M [16]: 



T 



and 



M R 0) = M ] , M { ° )+ = -H= h 22 , M| 0) = (7) 





2 Thcrcforc the lightest sterile neutrino in the vMSM is a candidate for the DM particle. 
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In this approximation we have two massive sterile neutrinos with equal mass M, the 
third neutrino is massless, and all active neutrinos have zero mass. It contradicts observable 
data [2, 3]. To adjust it next small terms are added to the matrix M DM [16]: 



m\\e ia m\2 m\3 



m£ ] = AM = ( m i2 m 22 e" i/3 

mi 3 m 33 e-^ 




h n 





hi 3 







h 23 


hi 





h 3 3 



M« + = ^ h 21 h 2 3 ,MP=0 (8) 



This correction violates U(l) e x U(l)^ x U(1) T symmetry, leads to the appearance of the 
mass of the third sterile neutrino and takes off the mass degeneracy for two heavier sterile 
neutrinos. It's also leads to the appearance of the extra small masses of the active neutrinos 
and nonzero mixing angles among them. 

In the terms of the introduced corrections Lagrangian (2) is 

C ad = -h aI L a N^ - MN C 2 N3 - —T^NjNj + h.c, (9) 

where Nj are right-handed neutrinos in the gauge basis. 

In order to find the masses of the active neutrino one has to make the diagonalization 
of the matrix M DM . The diagonalization undergoes in two steps. Firstly, M DM matrix 
is reduced to the block- diagonal form via the unitary transformation [17] in the "see-saw" 
approach: 

M block ^wv-wJ-^ T ^- 1M - ° W M r -° V ( 10 ) 



Mr J V M heavy , 

where 

W = ( 1 -jf E l jl Ee+ ) , e = M R H., D « 1 (11) 

and Mught = — (Mb) t (Mr)~ 1 Md, Mheavy = Mr are the mass matrix of the active and 
sterile neutrinos respectively. Now each block of the matrix M DM may be diagonalized 
independently by the matrix 

(12) 

The mass matrix of the active and sterile neutrinos is diagonalized by unitary transfor- 
mation Uu 2 y. 

UfMnghtUx = diag(mi, m 2 , m 3 ), U^M heavy U 2 = diag(M 1 , M 2 , M 3 ). (13) 
There is a standard parametrization [2] for Ui( 2 )' 

( 

U m = 



Cl 2 C 13 c 13 s 12 
■si 2 c 23 - ci 2 s 23 si 3 e l<5 c i2 c 23 - si 2 s 23 si 3 e J<5 
\ s i2 s 23 - c 12 c 2 3Si 3 e lS -c i2 s 23 - si 2 c 23 si 3 e l<5 






e ia ^ 2 |, (14) 
1 



1 



where Cy = cos^-, s^- = sin#y, 12 , #13, ^23 ar e the three mixing angles; 5 is the Dirac phase, 
and a%, a 2 are the Majorana phases. The angles can be in the region < 6^ < n/2, phases 
8, a^,a 2 vary from to 2tt. Each of the matrices U\ and U 2 contains its own, independent 
angles and phases. 

Then the elements of the Mu g ht can be defined by masses and elements of mixing matrix 
U of the active neutrinos: 



[M light ] a(S = m x U* al U; x + m 2 U* a2 U; 2 + m 3 U* a3 U; 3 . (15) 
The data that come from the neutrino oscillation experiments are presented in Tab.l: 



central value 99% confidence interval 
Am 2 21 = (7.58 ± 0.21) ■ 1(T 5 eV 2 (7.1 - 8.1) • 10" 5 el/ 2 
|Am| 3 | = (2.40 ± 0.15) • 10~ 3 eV 2 (2.1 - 2.8)10~ 3 eT/ 2 
tan 2 6 12 = 0.484 ± 0.048 31° < 12 < 39° 
sm 2 26 23 = 1.02 ± 0.04 37° < 9 23 < 53° 
* sin 2 20 13 = 0.11 (gig = 10°) 

Table 1. Experimental constraints on the parameters of active neutrinos [3], * — results of 
T2K Collaboration [13]: 0.03 < sin 2 29 13 < 0.28 in the case of the normal hierarchy and 0.04 < 
sin 2 2#i3 < 0.34 in the case of the inverted hierarchy. 



On the other hand, from the "see-saw" formula (in the approximation when the elements 
of the first column of the Yukawa matrix are neglected and M 3> m^) one can immediately 
obtain, that the mass of the lightest sterile neutrino is zero and the mass matrix of the active 
neutrinos has the form [16] 

v 2 

[M light ] a p = -^-^{h a2 hp 3 + h a3 hp 2 ), (16) 

and its eigenvalues is 

n v\F 2 F 3T \h + h\ 23 ) 
m a = 0, m (b ^ = — , (17) 



where F 2 = (h + h)n, m a is the mass of the lightest active neutrino, m c is the mass of the 
heaviest active neutrino. The sum over the neutrino masses is given by 

v 2 F 2 F 3 J-* , s 

2 s = V mi , 18 

i=i 

The system (16) has infinite number of solutions. Indeed, the replacement h a2 — > zh a2 , 
h a 3 h a3 /z (z is an arbitrary complex number) does not change the system. Then one can 
define the real quantity e 

e = F 3 /F 2 , e = \z\. (19) 
as an independent parameter of the model. 
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As it was shown in [18], the system (16) has good solutions for ratios of the elements of 
second column of the Yukawa matrix: 



.4 



M 



12 



12 



M- 



22 



.4 



M 



13 



13 



.4 



M 33 
M 23 



23 



M- 



1 ± 



l± Wl - 



1 ± a/1 - 



M n M 22 
M 12 2 

M n M 3 3 
M 13 2 



(20) 



33 



M 22 M; 



33 



M 23 2 



where A 12 = h 12 /h 22 , A n = hi 2 /h 32 , A 23 = h 22 /h 32 and M u is elements of Mn ght matrix. 
The ratios of the third column elements of the Yukawa matrix are expressed through the Aij 
elements: 

h 23 _ . M 22 h 33 _ M 33 

h^~ Al2 ^ i^-^w,- (21) 

Though formally there are eight different choices for the solutions (20), only four are 
independent. For example, if we fix the sign before the square roots in the expressions for 
A\2 and A\ 3 then A 23 is unambiguously determined by the relation 

A 23 = A l3 /A l2 . (22) 

The solutions (20) allow one to find the ratios of the elements of Yukawa matrix [18]: 



(hi2] h 22 ; h 32 ) 



Jarg(h 12 ) 



(h 13 ; h 23 ; h 33 ) 



y/l + \A 12 \ 

e iarg{h is ) 



-2 + |4. J -2 



(l; A l2 ; A 13 ) 



l 13| 



1 + \ A ^M^\ + 1^1317^1 



M 22 M 33 

a i2yt~ j A i3 



'M 



ii 



'Mi 



n 



where phases of h\ 2l h\ 3 are connected by condition 

arg(h 12 ) + arg(h 13 ) = arg(M xl ] 



(23) 
(24) 

(25) 



This is the exact solution of (16) that definitely expresses ratio of the elements of the Yukawa 
matrix via parameters of the active neutrino mass matrix. For fixed values of the active 
neutrino parameters there are only two choices for placing of the signs in the expressions 
for Ai 2 , Ai 3 , A 23 (20) which are not inconsistent with condition (22). These two variants are 
distinguished from each other by simultaneous replacement of the sign in front of square 
roots in the expressions for Ai 2 , A\ 3 , A 23 . It can be shown that such replacement of the signs 
leads to interchanging and conjugating of the ratios of elements of the second and the third 
columns of the Yukawa matrix, notably h 22 jh\ 2 -H- hl 3 /hl 3 , h 32 jh 12 -B- h\ z jh\ z [18]. 

As it was announced in Introduction, only the two heavier sterile neutrino take part in 
the production of the leptonic asymmetry. Therefore we will exclude the lightest sterile 
neutrino from consideration, so hereinafter indexes /, J take the value 2 or 3 referring to the 
two heavy sterile neutrinos. In this case there are 1 1 additional parameters in the vMSM as 
compared with SM. Seven of them we will identify with the elements of the active neutrino 
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mass matrix (m 2 , m 3 , 6 I 12 , #13, 6 23 , 8, ^2)- The other 4 we will define as follows: the average 
mass of two heavier sterile neutrinos M = M±Mi ; their mass splitting AM = M3 ~ M2 ; the 
parameter £ and the phase £ = arg{hi 2 ). 

Thus, we can parameterize the Lagrangian (9) in the following way: 



C 



1 - ~ y/e - ~ 
— — -=h a2 L a N a + -=-h a3 L a N 3 

J H 2\/ £ ^3 



- MN%N 3 - -AM ( N%N 2 + N£N 3 



h.c, (26) 



where a a j = h aI / Fj are defined by equations (23) and (24). 

Lagrangian (2) can be written in another basis, namely when the mass matrix of sterile 
right-handed neutrinos is diagonal. In this case the Lagrangian is 



■•ad 



-g aI L a N' I $- 1 ^N ,c I N , I + h.c., 



(27) 



where N'j are right-handed neutrinos and g a i are elements of the Yukawa matrix in this 
basis. 

Transition from presentation of Lagrangian (2) in gauge and mass basis can be made 
with unitary transformation that transfers mass matrix of right-handed neutrino to diagonal 
form [14, 16]: 



* I AM M 
M AM 



M- AM 
M + AM 



v = — (~ l 1 

V2\ 1 1 



So, the transition can be made by 

Nj = VjjN'j, g aI = h aJ Vjj. 



(28) 



(29) 



With help of this relations it will be useful to express Lagrangian (26) in terms of right- 
handed neutrino functions of Lagrangian (27) 



C 



ad 



2v 2 



% ^j=r - i\/Ea a3 j L a N' 2 + + \fea a3 ) L a N' 3 



— - (M — AM) N%N' 2 -~(M + AM)N' C 3 N' 3 . (30) 
After comparing (30) and (27) one can express Yukawa couplings in different presentations 



9oa 



9a3 



2v 2 



%a a2 

7^ 



i\fea a3 



da2 r- 

—f= + v ea a3 



(31) 
(32) 



The mass eigenstates neutrinos for Lagrangian with the mass matrix M DM (5) can be 
easily expressed through the states of neutrino of Lagrangian (27), particularly: 



N c 



-ee + N' c + eu L ~ N' 



(33) 



where N are mass eigenstates of the right-handed neutrinos in which they are produced and 
decay, v L are the active neutrinos of the SM in flavor basis, 

Sal - 8or = ^fg (34) 

is the mixing angle (e a j < 1). 

3 The computation of the leptonic asymmetry. 



As it was pointed in Section 1, leptonic asymmetry in the vMSM is generated due 
to decays of the heavier sterile neutrinos on SM particles. At temperature T -C T EW the 
interaction of the sterile neutrinos with SM particles via neutral Higgs field can be neglected. 
The only possible way of interaction of the sterile neutrino with matter is through the mixing 
with active neutrinos (33). 

For the sterile neutrino with the mass m, < Mj < 2 GeV the channels for the decay into 
two-body final state are: 

Ni ->■ n°u a ,n + e~,n~e^,K + e~,K~e^,r]u a ,r]'iy a ,p°u a ,p + e~,p~e^. (35) 

The channel of decay iV^ — > N\ + ... is strongly suppressed because of the small Yukawa 
coupling constants of N\. The decay of the sterile neutrino into the K° state is forbidden, 
because the composition of K° (ds) can not be obtained by decay of Z-boson. 

The three-body final state can be safely neglected and also the many hadron final state 
[10]. This last decay channels contribute for less than 10% for Mj < 2 GeV. For m n < Mj < 2 
GeV the decays into D-meson can also be neglected because its mass is not much smaller 
than 2 GeV. 

Let us consider the decay of the sterile neutrino in the vMSM. Sterile neutrino oscillates 
into active neutrino that decay into Z-bozon and active neutrino (or M^ ± -boson and charged 
lepton) in accordance with the SM. Z-boson (or H /± -boson) hereafter decays into quark- 
antiquark pair, see Fig.l. Since kinetic energy of this quarks are small enough the quark 
pair will form a bound state. Since Mi < 2 GeV <C Mz(w) we can use low energy Fermi 
theory and shrink the heavy boson propagator into an effective vertex and use for final state 
a meson, see Fig. 2. 

The process of sterile neutrino decay into charged lepton and charged meson through 
W /± -boson is described by charged current interaction 

£c = ^(j^)V CC , (36) 

where j„ = j l u cc + j„ cc is charged lepton and hadron current, 

ll CC = E ^7,(1 - 7 5 H, Ju CC = E V lmdmlu{l - l 5 )u n . (37) 

a n,m 

The indices m, n run over the quark generation, a = e, p,r and V is Kabbibo-Kobayashi- 
Maskawa (CKM) matrix. Similarly, the process of sterile neutrino decay into active neutrino 
and neutral meson through Z-boson is described by neutral current interaction 

£ N = V2G F (3? c ) + f NC , (38) 



S 




Figure 1: The decay of a sterile neutrino via Z-boson and H /+ -boson (the cross on line of a 
sterile neutrino means an oscillation of a sterile to an active neutrino). 



where = j l u NC + j^ NC is active neutrino and hadron current, 



J, 



NC 



r 



2v 



hNC 



H {4(1 - 7 5 ) - 2q f sin 2 6 W ) /, 



(39) 



where sum over / means sum over all quarks, t z - is the weak isospin of the quark, q$ - 
is the electric charge of quark in proton charge units, notably £3 = 1/2, q$ = +2/3 for u, c, t 
and t{ = —1/2, qj = —1/3 for d, s, b quarks. 

The matrix element corresponding to Feynman diagram of sterile neutrino decay (see, 
e.g., Fig. 1,2) can be obtained from the interactive effective Lagrangian [14]. For example, 
effective Lagrangian of decay of I sterile neutrino into the 7r ± , 7r° final states is: 



^eff 



G } 



■M I Ue aI u a {l+^)N I ^ 



+ 



^MjUV ud e aI e a ((1 + 75 ) - ^(1 - 75) ) Njn~ + h.c. 



(40) 



where Gp is Fermi coupling constant, Mj is the mass of J-sterile neutrino, m a is the mass 
of the charged lepton of a generation, /„. is the 7r-meson decay constant that is defined as 



(n+\u(l+^) lu d\0) = -f n -( Pn )„ 
where p n is the pion 4-momentum. 



(41) 




Figure 2: Effective low-energy decay of a sterile neutrino into 7r° meson and active neutrino. 
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a) 




N i 




7t, 71 



b) 



A/, c 




Figure 3: Example of one-loop diagrams of the decay Nj — > z/ Q 7r c 



The leptonic asymmetry e can be defined as 



N- 



+ r 



(42) 



where T^^i is the total decay rate of sterile neutrinos into leptons and T N ^i is the total 
decay rate of sterile neutrinos into antileptons. 

At tree level the decay rates of the sterile neutrinos into leptons and antileptons are 
equal. Therefore we must compute the one loop diagrams, see Fig. 3. In the case of nearly 
degenerated sterile neutrinos the contribution from the diagrams presented at Fig.3&) can 
be neglected as compared with diagrams presented at Fig.3o). Indeed the propagator of 
the sterile neutrino in the diagrams a) type is proportional to 1/AM in the center of mass 
frame. The leading order contribution to the leptonic asymmetry comes from interference 
between one- loop diagrams and tree-level diagrams [19]. In this case T N ^i — T N ^f ~ O 4 and 
T N ^[ + T N ^j ~ O 2 , the leptonic asymmetry is suppressed. 

In our case, when the mass splitting between the two heavier sterile neutrinos is very small 
and it is of the same order as their decay rate (we obviously will see it later), the oscillations 
between Nj and Nj are important, see Fig. 4. So, the corresponding mass eigenfunctions 
are no longer the Nj states, but a mixture of them, namely ipj [20, 14]. It is these physical 
eigenstates which evolve in time with a definite frequency. The subsequent decay of these 
fields will produce the desired lepton asymmetry 



r - ' 



(43) 



where and are the total decay rates of the sterile neutrino mass eigenfunctions 

ipi into leptons and antileptons correspondingly. In this case the leading order contribution 
to the leptonic asymmetry comes from tree-level diagrams. 

In general case the correct description of the processes can be made in frame of the density 
matrix formalism, see, e.g., [5]. We will follow a simpler way by considering a non-hermitian 
Hamiltonian. The effective Hamiltonian in the basis of {^2, N 3 } is the H = H + AH, where 
Hq is the diagonal Hamiltonian of equal mass particle 



H 



M 
M 



(44) 
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71, 71 



v„, e 




Nj 



7t, 71 



Figure 4: Contributions to the effective Hamiltonian. 



The corrections to this Hamiltonian are given by the one-loop diagrams, see Fig.4: 



AH 



-am - |r 2 

i -p 



AM 



23 

■ir 3 



(45) 



The dispersive part of these diagrams can be absorbed in the mass renormalization of the 
fields [20] and it brings to appearance of the mass splitting AM. The absorptive part of the 
diagrams will define total decay rates of the sterile neutrino Tj and the rate of oscillation 
between sterile neutrinos T 2 3. 

Total rates of /-sterile neutrino decays into charged mesons and leptons of a-generation 

are 



_ Glft\V ud \ 2 M* 
8tt 



Q aI \ 2 S{M, m a ,m v ) 



M 2 



M 2 



1 + 



M 2 



(46) 



1 O.K 



11) 



G 2 F f K \V us \ 2 M s 



Stt 



Q aI \ 2 S(M, m a ,m K ] 



M 2 



in 



A 



M 2 



1 + 



M 2 



(47) 



r; 



GW P \V ud \ 2 M* 



4irm 2 



e aI \ 2 S(M,m 



M 2 



m P 
M 2 



mz 



2m 2 / 



M 2 



(4* 
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where 

/ / {m — m n ) 2 \ ( (m + m„) 2 \ , . 

S(M„m a ,m) = '—Jff^) (l - StT^J' (49) 

and values of decay constants and elements of CKM matrix are given in [2]: /„- = 0.131 GeV, 
f K = 0.16 GeV, g p = 0.102 GeV 2 , \V ud \ = 0.97, \V US \ = 0.23. 

Total rates of /-sterile neutrino decays into neutral mesons and active neutrinos are 

rr° = m -» *° + „j = (i - g) 2 , (50) 

r- = iw -> Po + „j = -^-|e^|« (i + 2^) (1 - ^) , (51) 
r- = r(jv, - n + *J = ^jf^l 2 (1 - . W 

Tf = T(N, -» r/ + „J = (l - g)\ (53) 

where /„ = 0.156 GeV, f v = -0.058 GeV [2]. 

As one can see the decay rates into p^p mesons are slightly different because they are 
vector mesons. The adduced decay rates (50) - (53) were obtained in [21, 10]. The total 
decay rate of sterile neutrino decay into mesons and leptons is sum of the rates over all decay 
channels A (35) and over leptonic generation: 

r / = J]r J oA e(y aA ) ) (54) 

a, A 

where y\ is the difference of the I sterile neutrino mass and total mass of all final particles 
of the decay channel A; Q(x) is the usual Heaviside function. The rate of oscillation between 
/ and J sterile neutrinos (Tjj) can be expressed through the decay rates 

a, A 

The eigenvalues and corresponding eigenfunctions of the non-hermitian Hamiltonian H = 
H + AH are given by 

U2 = M- l -(T 2 + T 3 )-^c, 

w 3 = M-|(r 2 + r 3 ) + - c , 

where N is a normalization factor and 

c = v / (4AM-z(r 3 -r 2 ) 2 -4(r 23 ) 2 , b = (4iAm + (r 3 - r 2 ) + zc)/(2r 23 ). 

It should be noted the sterile neutrinos are not initially in the state ip2 and ips, but in 
the state iV 2 and iV 3 . The fact is that sterile neutrino where in thermal equilibrium before 
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they propagated freely. The equilibrium was maintained by the weak interaction between 
the sterile neutrinos and particles in the background. The weak interaction eigenstates are 
N 2 and N 3 , therefore at the beginning the sterile neutrinos are in the state N 2 or N 3 . In 
general the initial state of sterile neutrino is the superposition of N 2 and N3 states and can 
be described by a density matrix: 



Pinitiai = p(t = 0)=J2 aj|Wj(0))(jV/(0)|, (5* 



1=2,3 

where a 2 + 0:3 = 1. It was shown in [1 1] that leptonic asymmetry dependence on parameter 
aj can be neglected. We confirmed this statement and, hereafter, we will consider the 
symmetric initial state a 2 = «3 = 1/2. 

The time evolution of the density matrix can be obtain in a simple way. Since 

\^i) = Uu\Nj), (59) 

where 



'N 

the time evolution of \Ni) state is known 

Wi(t)) = C/^e-^I^O)) = U^e-^UKjlNjiO)) = Ru\Nj(0)). (61) 

Thus 

3 3 



P® = \ E Rm{tYRu{t)\NM){N K (S))\ = \ Y^{R ] R)kj\NM)(N k {Q)\. (62) 

I,J,K=2 J,K=2 

The average production rate of leptons is given by 



/*oo p 1 poo i* 

r =J dtjdU 2 J2Tr[\l)(l\p(t)} = -JdtjdU 2 Tr 



J2(R f R)Kj(l\Nj(0))(N K (0)\l) 

\,K,J 

-\ fdtfdU^i^R^KjAjtA^, (63) 

ZJo J LJ,K 



where sum over I means sum over all leptons generations and include charged leptons and 
active neutrinos, (l\Nj(0)) = Aji is the transition amplitude of the decay of / sterile neutrino 
into a lepton at tree level that includes all possible channels of reaction, and dU 2 is the 
differential 2-body phase space 

where p, q, k are 4-momentums of initial and final particles in decay. 
Similarly the production rate of antileptons is 

f = \ fit [dU 2 {R ] R)kjA* ji A kj . (64) 

ZJo J LJ,K 
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The measure of the leptonic asymmetry is given by 
A r-T fdtf dU 2 Im((R^R) 32 )Im(A* 2l A 3l ) 



The integration over dU 2 gives [14]: 

/ dtIm((tfR) 32 ) J2 a Im(Q* a2 Q a3 )V a 



I ^Ea((^)2 2 |e a2 | 2 + (WR) 33 \e a3 \ 2 + 2Re(e* a2 e a3 )Re(Rm) 23 )v a 



(66) 



where V a is defined via sum over all possible channels of sterile neutrino decays into leptons 
of generation a 

"p aA 

4 The restrictions on the parameters of the vMSM 



As it was pointed in Section 1, the leptonic asymmetry of the Universe has to be con- 
strained by condition (1) at the moment of the beginning of the DM particles production. 
It allows us to constrain parameters of the vMSM . To do it, we can construct the leptonic 
asymmetry (66) as function of only three parameters of vMSM: M, AM, e. 

We do it in the following way. Leptonic asymmetry function (66) is maximized over 
phases 5, a 2 , £ (and a\ in case of the inverted hierarchy) and is taken at central value of 
active neutrino mass matrix parameters 3 , see Tab.l. This function contains dependence on 
ratios of the Yukawa matrix elements in mixing angle Q / (34) that can be expressed through 
solutions (20) with two possible choice of sign consistent with condition (22). So far as the 
relation for leptonic asymmetry (66) has no symmetry for interchanging and conjugating of 
the ratios of elements of the second and the third columns of the Yukawa matrix we have to 
consider two variants of the solutions. For fixed values of the mixing angles and phases we 
will designate allowed solution of (20) with 2 or more sign (+) as solution of A type, and, 
vice-versa, the solution with 2 or more sign (— ) we will designate as solution of B type. It 
should be noted that our results (23), (24) for B type of solution coincide with results of [8] 
where the ratios of the elements were obtained in the particular case 8 r3 — > 0, 9 23 — > 7r/4. 
We separately consider the case of 6 13 = and 9 13 = 10° also. 

Thereby we construct allowed regions (A > 10 -3 ) in plane of parameters AM and e at 
fixed values of M. 

For the case of the normal hierarchy the deference between the case of #13 = or # 13 = 10° 
and the case of solution of A or B types is not essential, so we illustrate allowed regions with 
help of only one figure on Fig. 5. For the case of the inverted hierarchy, the difference between 
the case of solution of A or B types is not essential, but the cases of 6*13 = and 6*13 = 10° 
are substantially different. So we illustrate allowed regions with help of two figures on Fig. 6. 

It should be noted that we investigated form of the allowed regions not only for the central 
value of #13 angle, but for range given by data of [13]. We conclude that in case of the normal 



3 In case of the normal hierarchy we have mi = 0, m-i = \/ Am^ = 0.009 eV, = yj\ Am^l + Am| 



2.1. 



0.05 eV. In case of inverted hierarchy we have mi = -^Am^l — Ato^i — 0.048 eV, m-i = \J\ Am 
0.049 eV, m 3 = 0. 
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3.0 -2.5 -2.0 -1.5 -1.0 -0.5 



Figure 5: The grey areas are the regions of parameters where A > 10 3 for the case of the 
normal hierarchy. The areas correspond to M — 0.3 GeV (bottom), M — 1 GeV (middle), 
M = 2 GeV (top). 

hierarchy the regions are almost not sensitive to value of 6*13 in range < #13 < 16°. In case 
of the inverted hierarchy it is true for the regions on Fig6 b) and 6*13 < 18°, but for #13 = 
the allowed regions are appreciably different. 

Also we illustrate regions where maximum of A can be more then 2/11 on Fig. 7 (white 
inner figures) for the case of both hierarchies. We do it only for the mass M — 1 GeV 
because this regions are at small values of e and it will not intersect with other subsequent 
constrains. Moreover, at some values of phases leptonic asymmetry in this region can be less 
then 2/11 and so we can not exclude this region ultimately. 

By way of example, we present possible values of / sterile neutrino decay rate Tj (54) and 
rate of oscillations between / and J sterile neutrinos Tjj (55) for M — 1 GeV and #13 = 10° 
on Fig. 8. As one can see the values of T/, T/j are really of the same order as AM. It confirms 
previous assumption about necessity of taking into consideration oscillations between sterile 
neutrinos. 

In order to create a leptonic asymmetry the sterile neutrinos should be out of thermal 
equilibrium. That means that 

T 2 < H, (68) 
where H is Hubble parameter that determines the expansion rate of the Universe. In the 




Figure 6: The grey areas are the regions of parameters where A > 10 -3 for the case of the 
inverted hierarchy. The areas correspond to M = 0.3 GeV (bottom), M — 1 GeV (middle), 
M = 2 GeV (top). Figures a) and b) represent the case of 6* 13 = 0° and 10° correspondingly. 
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Figure 7: The grey areas represent regions of parameters where 10 3 <A<2/11 for M — 1 
GeV in case of normal (a) and inverted hierarchy (b). 

radiative dominated epoch Hubble parameter is given by 

T 2 



H 



M 



(69) 



PL 



where Ml 



90 



1 pl — y 8-K 3 g*(T) MpL, Mpl = 1.22-10 19 GeV is the Planck mass, g*(T) is the internal 

degrees of freedom [22]. At temperature T ~ 1 GeV we can take g* ~ 65. 
So we get condition 

VWl^ < T. (70) 

The out-of-equilibrium condition means that sterile neutrinos should decay at a temper- 
ature smaller than their mass (T < M). Moreover the sterile neutrinos should decay before 
the creation of DM so that the leptonic asymmetry enhances the DM production. The DM 
is created at T ~ 0.1 GeV. Therefore, 



n I < X £ < 

' ~ lGeV ~ 



M 



lGeV 



(71) 



We illustrate on Fig. 9 the region of values of M and e where condition (71) is satisfied 
for the case of the normal hierarchy. At scale of parameters presented on Fig. 9 the difference 
between the case of #13 = or 9 13 = 10° and between the case of solution of A or B types is 




r- -20 



-22 



-24 



-26 




-21 



-23 



-25 



Figure 8: The values of rates Lo g 10 (T 1 / lGeV) , Logio(r 2 3/lGeV) and Logio(AM/lGeV) for 
leptonic asymmetry A > 10 -3 for the case of M = lGeV and #13 = 10° are on the ordinate 
axis: a) the case of normal hierarchy (A type of solution), b) the case of inverted hierarchy 
(B type of solution). 
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Figure 9: The case of the normal hierarchy. The points on grey and red regions satisfy 
constraint (71). The region on the right from vertical red line satisfies condition (1) also. 



small, so we present only one figure. It is not true for the case of the inverted hierarchy, see 
Fig.10. 

It should be noted that region on Fig. 9 is almost not sensitive to value of #13 at range 
< #13 < 16°. In case of the inverted hierarchy it is true for the regions on Fig.10 b) and 
0i3 < 18°, but for the case of #13 = the regions are appreciably different. 

As one can see there are regions on Fig. 9 (red) and Fig.10 a) (red and blue) where 
conditions (1) and (71) are satisfied simultaneously. This region of parameters is suitable for 
DM production in the vMSM . For the case of inverted hierarchy and nonzero value of #13 
we have no region that is suitable for DM production. So, in v MSM for physical nonzero 
value of #13 and mass of sterile neutrino m n < M < 2 GeV DM production can be realized 
only in case of normal hierarchy of active neutrino mass. 

The region suitable for DM production (the case of the normal hierarchy and nonzero 
#13) can be used to obtain constraints for mass splitting of the sterile neutrino. Fixing mass 
of the sterile neutrino one obtains possible values of e (see Fig. 9) and using Fig. 5 one can 
obtain possible values of the mass splitting for the sterile neutrino with mass M. If mass 
of the sterile neutrino is on lower boundary of the allowed mass range (M ~ 1.4 GeV) than 
value of AM is exactly known (AM « 5 • 10~ 21 GeV). If mass of the sterile neutrino is on 
upper bound of the allowed mass range (M = 2 GeV) than AM can possess the values from 
the range 10" 21 < AM/lGeV < 10" 20 . 

Some existing experimental data restrict the area of parameters of vMSM. For M < 0.45 
GeV the best constraints come from the CERN PS191 experiment. For 0.45 < M < 2 GeV 
the constraints come from the NuTeV, CHARM and BEBC experiments. The range of 
parameters admitted by these experimental data is summarized in [23]. These parameters 
are the mixing angle (0 + O)22 (it defines the range of reactions with sterile neutrino) and 
the mass of the heavier sterile neutrino M. 

To compare obtained in the present paper constraints on the vMSM parameters (see 
Fig. 9 and Fig.10) with constraints summarized in [23] one has to rebuild allowed regions in 
the space of parameters M and 9„ N = (0 + 0) 2 2- 

In general case the relation between (0 + 0)22 and e is quite difficult. Really, in accordance 
with (29) and(34) we have 



(0+0) 




(V + h + hV) 




(F 2 2 + F 2 -2|/ i + / i | 23 sinx), 



(72) 
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Figure 10: The case of the inverted hierarchy: a) #13 = 0, b) #13 = 10°. The pink, grey and 
red regions corresponds to the A type of solutions. The grey, red, sky blue and blue regions 
corresponds to the B type of solutions. The points on this regions satisfy constraint (71). 
The region on the right from red line satisfies condition (1) also. 



where x — ar 9[(h + h) 23]- Using (17) - (19) we get 



M 
u 2 e 



:m. 



m b 



eFo 



\h + h\ 23 = —^{m c 



1^ 



m b 



and 



(0+0) 



22 



2Ms 



2 m c -m b 

e —2 esmx 

m c + m b 



(73) 



(74) 



The problem is in parameter \ that is a complicated function of many parameters. But for 
our case (see Fig.9 and Fig. 10, e < 0.16) we can use approximate relation 



(0+0) 



22 



m c + m b 
2Mb 



(75) 



The imposition of our constraints are presented on Fig.9 and Fig. 10 for nonzero value 
of #13 and summarized constraints from [23] is presented on Fig. 11. Above the line marked 
"BAU", baryogenesis is not possible: here sterile neutrinos come to thermal equilibrium 
above the Tew temperature. Below the line marked "See-saw", the data on neutrino masses 
and mixing cannot be explained using "see-saw" mechanism. The region noted as "BBN" is 
disfavoured by the considerations of Big Bang Nucleosynthesis. The region marked "Experi- 
ment" shows the part of the parameter space excluded by direct searches for singlet fermions. 
The regions market "Cos" , "A" and "DM" were builded in this paper. The grey and blue 
region "Cos" shows the parametric space allowed by cosmological constraint (71) (grey re- 
gion corresponds to A and B type of solution, blue region corresponds to B type of solution), 
the dashed region market "A" shows the parametric space allowed by constraint (1), the 
red region marked "DM" shows the parametric space where constraints (1) and (71) are 
noncontradictory. The last region is preferred for DM production according to calculations 
of the present paper. 

The red region marked "DM" is shown on Fig. 12 in the scaled- up form. The difference 
between the case of #13 = 0° or #13 = 10°, and between type of A or B solutions is illustrated. 
As one see the choice of solutions of A or B type makes greater change in the allowed region 
then the choice of #13 = 0° or #13 = 10°. 
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Figure 11: The imposition of our constraints and summarized constraints from [23]: a) the 
case of the normal hierarchy, b) the case of the inverted hierarchy. 

5 Conclusion 

In the present paper we consider the leptonic asymmetry generation at T <C Tew when the 
masses of two heavier sterile neutrinos are between m v and 2 GeV. 

We conclude that oscillations and decays of sterile neutrinos can produce a leptonic 
asymmetry that is large enough to enhance the DM production sufficiently to explain the 
observed DM in the Universe, but only for the case of the normal hierarchy of the active 
neutrino mass. The allowed range of parameters is narrow and it is presented on Fig. 11 and 
Fig.12. It should be noted that allowed mass range for heavier sterile neutrino is 1.42(1.55) < 
M < 2 GeV for B (A) type of solutions and the mixing angle between active and sterile 
neutrino is -7.91(-7.98) < Log lo (e + 0) 22 < -8.41(-8.35) for B (A) type of solutions. If 
mass of the sterile neutrino is on lower boundary of the allowed mass range than value 
of AM is exactly known (AM m 5 ■ 10~ 21 GeV). If mass of the sterile neutrino is on 
upper bound of the allowed mass range than AM can possess the values from the range 
10~ 21 AM/lGeV < 10 -20 . For the case of the inverted hierarchy there is no region 
suitable for DM production. 

The big range of parameters of the vMSM is not forbidden by the existing experimental 
data, see Fig. 11. Combining of this range with our constraints (red region "DM" on Fig. 11) 
leads to conclusion that improvement of previous experiments, as NuTeV or CHARM, of 
one or two order of magnitude can exclude the vMSM with M < 2 GeV or detect the 
right-handed neutrinos. 

It should be noted that our constraints are quite a rough and can be used only for 
estimation. Really, the form of red region " DM" is very sensitive to cosmological constraints. 
Applied condition O.lGeV < \jMpi,*Ti < M is very approximate. The correct description of 
the processes can be made in frame of the density matrix formalism or Boltzmann equations. 
Our computation is not valid for M > 2 GeV. However, the extrapolation of our result, see 
Fig.ll, suggests that the range of admitted parameters for the case of the normal hierarchy 
becomes bigger for masses above 2 GeV. We expect that for masses above 2 GeV DM 
production can be realized for the case of inverted hierarchy too. 

During computation we used two types (A or B) of solutions (20). This is due to the fact 
that ratios of the Yukawa matrix elements (enter into the expression for the mixing angle 
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Figure 12: The red region "DM" from Fig. 10 in the scaled-up form: a) 6* 13 = 0° type A 
(dark) and type B (light); b) type A: 13 = 0° (white) and 6 13 = 10° (black); c) type B: 
6*i3 = 0° (white) and #13 = 10° (black). The variable M/l GeV is along the abscissa axis 
and the variable Log 10 (9 + 6)22 is along the ordinate axis. 

Q a i) can be expressed through solutions (20) with two possible choice of sign consistent with 
condition (22). It is closely related to the symmetry of (16) under replacing the elements of 
the second column of the Yukawa matrix by elements of the third column. This two variants 
are equal in rights. 

The computation of the leptonic asymmetry in the applied simple model allows us to make 
some conclusions that, seemingly, will be correct and under more rigorous consideration. 
Namely, the initial state of the right-handed neutrino in form (58) are not important for 
lepton asymmetry generation (the final results are not sensitive to values of the constants 
ai). For the case of normal hierarchy the deviation of the mixing angle 6*13 from its zero 
value (up to value 16°) almost does not change the region suitable for DM production. For 
the case of inverted hierarchy results are different for #13 = and #13 7^ 0. Our calculations 
indicates that case of 6* 13 = leads to existing of region suitable for DM production, but at 
nonzero values of 6^3 this region does not exist. Values of 6*13 in range 6> 13 < 18° (# 13 7^ 0) 
almost does not change the region suitable for DM production. 

It's essential to note that during computations we have used functions maximized over 
unknown parameters of the model (phases <5, CK2, oc\) . If the maximization procedure was 
not performed the final functions are sensitive to values of mentioned phases. So, the ob- 
tained results are very optimistic. But if the proposed on Fig. 11 region of parameters "DM" 
will be forbidden by experiment data it will mean that mass of heavier sterile neutrinos must 
be lager 2 Gev. 

An essential assumption we have made is that the background effects are negligible. 
We do not have justify that it can be neglected in the thermal bath of the universe. For 
simplicity the computations were made at zero temperature. A rigorous justification of this 
assumptions is needed. 

It should be noted that region suitable for DM production in vMSM was recently cal- 
culated in frame of more general formalism in [24]. Certainly, results of [2 1] somewhat differ 
from our simple calculations. 
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